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We develop a quantum noise approach to study quantum transport through nanostructures. The nanostruc- 
tures, such as quantum dots, are regarded as artificial atoms, subject to quasi-equilibrium fermionic reservoirs of 
electrons in biased leads. Noise operators characterizing the quantum fluctuation in the reservoirs are related to 
the damping and fluctuation of the artificial atoms through the quantum Langevin equation. The average current 
and current noise are derived in terms of the reservoir noise correlations. In the white-noise limit, we show that 
the current and current noise can be exactly calculated by the quantum noise approach, even in the presence of 
interaction such as Coulomb blockade. As a typical application, the average current and current noise through a 
single quantum dot are studied. 
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I. INTRODUCTION 

Quantum transport through nanostructures is of importance 
in nano-science and nano-technology. Many electronic de- 
vices based on nanostructures, such as single-electron tran- 
sistors, have been studied in the past decades for their poten- 
tials in various applications. Recently, in the efforts aiming at 
the coherent control of single electrons or electron spins, the 
quantum transport methods have been used to detect the quan- 
tized motion of electrons in nanostructures.^ Besides the av- 
erage current, the current noises also contain useful informa- 
tion about the quantum dynamics in nanostructuresi^ii^ 

Various theoretical approaches have been developed to treat 
the quantum transport problem. The Landauer-Biitikker for- 
mula has established the basic relationship between scattering 
amplitudes and currents through nanostructures.^ The non- 
equilibrium Green's function method provides a perturba- 
tion scheme to deal with the many-body interaction effects in 
quantum transport^ In the past few years, approaches based 
on notions in quantum optics were developed to study time- 
dependent quantum transport processes in solid-state struc- 
tures. Most of these quantum optics approaches adopt the 
density matrix formalism, with master equations or rate equa- 
tions in the Schrodinger picturei 10 ' 11,12 ' 13,14 ' 15 ' 16 ' 17 ' 18 Very re- 
cently, quantum Langevin equation in the Heisenberg picture 
was also used to establish the quantum rate equations to study 
the transport problem^ 

In this paper, instead of the master equations or the rate 
equations, we will develop a quantum noise approach based 
on the quantum Langevin equation to the quantum transport 
problem. Essentially, we recognize that a general quantum 
transport problem can be regarded as a system-plus-reservoir 
problem. In this sense, the total system is divided into several 



sub-systems (see Fig[T). The central system {system for short) 
is a nanostructure, such as a quantum dot or coupled quantum 
dots. This subsystem contains several discrete electronic en- 
ergy levels, resembling an artificial atom. The electrons in the 
leads, which have a continuous energy spectrum and are kept 
in quasi-equilibrium, constitute the fermionic reservoirs. The 
electrons in the reservoirs can be treated as free quasi-particles 
with the screened Coulomb interaction taken into account as a 
renormalization of electron effective mass. The central system 
and the reservoirs are coupled together to each other through 
hopping across the barriers. With this observation, it is nat- 
ural to treat the quantum transport problem in the framework 
of the quantum open system method, the quantum Langevin 
equation, a standard approach in quantum optics to study cav- 
ity photon decay and atom damping. 

As compared to the application in quantum optics, the 
quantum Langevin approach in the quantum transport prob- 
lem has two features to be singled out: (i) The reservoirs con- 
sist of electrons, which are fermions while the baths in quan- 
tum optics are bosonic, and (ii) when finite biases are applied 
between different leads, the electronic reservoirs in different 
leads are in quasi-equilibrium with different chemical poten- 
tials but do not stay in equilibrium with each other. Our inves- 
tigation in this paper will refine these features. As illustrative 
applications of our approach, the resonant transport through 
a single quantum dot is investigated for both the single-level 
case and the Coulomb blockade case. 

The quantum Langevin approach is a natural formalism 
to study the noise spectroscopy of quantum dynamics in 
nanostructures, 5 which is particularly interesting for small 
quantum systems where the signals are often much weaker 
than the shot noises. When the coupling between the leads 
and the nanostructures can be described in the Markovian ap- 
proximation, which is justified in large bias cases, the quan- 
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FIG. 1: (color online) Schematic illustration of transport through 
a quantum dot. The whole system is divided into three parts, the 
central system, and the left and right reservoirs. The central system 
is characterized by eigen-states \i) with discrete energies. Quantum 
noise operators £.{t) and ■??(?) are introduced to describe the reser- 
voirs. 



turn Langevin approach provides an exact treatment of the 
interaction within the nanostructure. Furthermore, the quan- 
tum Langevin equation establishes a fundamental relationship 
and analogy between photon emission and electron tunneling 
processes, providing new understanding of quantum transport 
phenomena with notions and methods from quantum optics. 

The paper is organized as follow. In Sec. [TT] we introduce 
the basic concepts and the general formalism of the quan- 
tum noise approach to treat the quantum transport problem. 
In Sec. HlJand llVl we apply the quantum noise approach to 
transport through a single quantum dot containing a single 
level and double energy levels, respectively. In Sec. |VJ we 
show the relations between our approach and other quantum 
transport theories. We conclude and give an outlook of our 
approach in Sec lVIl 



II. GENERAL FORMALISM 

A. Quantum Langevin equations for quantum transport 

In general, the quantum transport problem of nanostruc- 
tures can be modeled by the following Hamiltonian: 

H = H sys ( ai ,al)+H lead + H T , (1) 

where // sys describes the nanostructure, such as a quantum 
dot, with multiple discrete energy levels. The leads, which 
play the role of reservoirs, are described by the Hamiltonian 
//lead- The electron tunneling between the leads and the nanos- 
tructure is included in Hj. For the two-lead case, the leads 
Hamiltonian Z/i ea d and the tunneling Hamiltonian Hj can be 



written as 

ttiead = J^tUo^blh + J^najf^Cj, (2a) 

* J 

H T = ih^^ ik b' k ai + ih^ZijCjCii + h.c, (2b) 

i,k i,j 

where b k and cj are the annihilation operators of the left and 
right leads with continuous spectra hal^ and %coy, respec- 
tively. The tunneling is characterized by the coefficients ft. 
and dj. Note that we have neglected the interaction in the 
leads as a common approximation for a Fermi sea with the 
Coulomb interaction effectively taken into the renormalized 
quasi-particle spectra. 

Now, we consider the Heisenberg equations of motion of 
the system and reservoir operators. For simplicity, we show 
equations of motion for the simplest single level case, i.e. 
a,(?) = a(t), ftt = ft, — £/and Z/ sys = fiuoa^a. The multi- 
level case will be discussed later in this paper. Straightforward 
calculation gives 



a (f) = -ico a - ^ ftfc* - ^ £jCj, (3a) 
* j 

b k (t) = -mfb k + ^ k a, (3b) 

cj(t) = -icj^Cj + ^ja. (3c) 



In the following, we try to eliminate the lead variables from 
the equation of motion (l3aT l of the system operator. To this 
end, the formal solution for b k (f) is written as 

b k (t) = (0) + ft J dt' [e-^'-^a {t')\ . (4) 

With this formal solution, the following relation is obtained 

2ftZ>*« = ~£in(*) + ^a(f), (5) 
k z 

where the input noise operator Xi n (0 due to the left lead is 
defined as 

£m(t) = -^ k e- iJ ^b k (0). (6) 
k 

The damping term in Eq. (O arises from the Markovian 
approximation? 1 ' under the continuous limit 

k 

= j dcj k [d (u> k ) f (co k ) e -H°('-<'> 

= y L s(t-f), (7) 

where D{<jL> k ) is the density of the states of the lead. Here we 
have assumed that D(<y^)^ 2 (wi) is flat around the frequency 
cjq and 

7L = 2nD(u )f(cJo), (8) 
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is widely used as the tunneling rate in nanostructure quantum 
transport problems. 

Similarly, for the right lead 



2?A-(t) = -% 1 (t) + ^fl(0, 



(9) 



where the noise operator of the right lead % (f) is defined as 

%n{t) = -^ j e- Ud f t c j (0), (10) 

and y R is the tunneling rate to the right lead. Using Eqs. (0 
and (0, we obtain the quantum Langevin equation for the sys- 
tem operator a{t): 

a (t) = -KUoa(f) - ^p^a (t) + £ in (t) + <R in (f) . (11) 

Similar to the cases in quantum optics, the two electronic 
leads, which play the role of fermionic reservoirs, induce the 
damping and the fluctuations through the noise operators. 

B. Projection operator formalism for interacting systems 

For an interacting system, the complexity of the quantum 
Langevin equations arises from the evolution induced by the 
system Hamiltonian // sys . To deal with such complexity, we 
introduce the projection operators of the interacting system in 
this subsection. 

Though there may be interaction between the electrons in 
the system Hamiltonian // sys , the artificial atom can always be 
considered as consisting of a few discrete many-body energy 
levels. In other words, we can diagonalize the system Hamil- 
tonian // sys as 

N N 

H sys (at, at) = ^ ha>i = ^ ftw,-o- s , (12) 

i=l ;=1 

where cr,y = is the projection operator from \ j) to |z), with 
|/) being the eigen-state of H sys of energy hu,. 

In general cases, determined by the system Hamiltonian 
i/ S ys, the fermionic operators a, and aj can be written in terms 
of the projection operators as 

at = J] 7*V« and a] = J] T®(r u , ( 1 3) 



where T ( '^ and T (,) me N x N matrices associated with the 
fermionic operators a, and at, and T {l) = (r (,) ) . 

The commutative relation between the projection operators 
CTij can be calculated with cr^o"a = (TuS^. In the following 
calculations we also need to identify the commutative relation 
between the projection operators cry and the reservoir opera- 
tors, i.e. bk{b' k ) and c y (cj). Note that the eigen-state \i) is also 
an eigen-state of the electron number operator in the quantum 
dot N = Yjidjcii, i-e. N\i) = N, Thus, the projection op- 
erator cry corresponds to a definite electron number change 



Ni — Nj, which is either odd or even. Consequently, cry and 
bk{b\) have the following (anti-)commutative relation: 



[o-ij, b^ = a-ijbk + gijb k o- u = 0, 
where the factor 



1, for Ni-Nj = odd 
1-1, f 'or Ni-Nj = even 



The Heisenberg equation for the cry is 

o-yW = ^K ys ,o-y] + -[fi r r,o-y] 

= -iAijCTjj(t) 

a.k.m.n 



(14) 



(15) 



- J] [(^<r„-h.c.),<r i J. (16) 

a,j,m,n 

With the help of the definition of noise operators and the first 
Markovian approximation, we obtain the quantum Langevin 
equation for the projection operator cry 

& u (t) = -;'A, 7 cry(f) 



m,m' 

+ J] C^Xin (t) CT mm , (f) + Uin -» %n) , (17) 



where Ay = o)j—a>i, and the coefficients DjjL,, CjL,, and C 
are defined as follows 



mm' 



D'L = A l i ml +Al n , + B'i m , + B , l J nml , (18a) 

cL = EffiW (18c) 



with 



mm' 


- \^ a- .f^fte) 
/ j & l J mi jm' ' 


(19a) 


A ij 

mm' 




(19b) 


mm' 


- A V T^T^ 

— "m'j 1 mv 1 vi ' 

v,a,a' 


(19c) 


mm' 


- x V 7'(o')'f(a) 
Y^a 1 


(19d) 



In principle, the quantum Langevin equation for the system 
operators is equivalent to a quantum stochastic equation if we 
introduce the quantum Wiener processes and the properties 
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of the their solution can be discussed by defining the quan- 
tum stochastic integration^ Thus, we point out that the quan- 
tum transport problem provides an experimentally accessible 
proving ground for the quantum stochastic theory. Instead of 
further discussing the mathematical properties of Eq. $1% , in 
this paper we will focus, through concrete models, on how to 
derive the observable quantities in quantum transport. 



C. Boundary Relation and Causality 

Besides the input noise operators Xin (?) and % n (f), the out- 
put noise operators^2 can be defined as 

Xout(f) = -Yj^e-^^hik), (20a) 

flout (0 = -^je-^'-'Icjitf), (20b) 
i 

where tf is a time in the remote future. Similar to Eqs. (O 
and ©, the first Markovian approximation gives the following 
relations 

J] {t) = -Xout (f) - (f) , (21) 

k L 

According to Eqs. (0 and ( |2TI ). the "boundary relation" be- 
tween the noise operators and the system operator is 2 ^ 

Xta (t) - Xout (t) = y L a(t). (22) 

Similarly, for the right lead 

%n(t)-Kut(t)=yRa(t). (23) 

According to the quantum Langevin equation ( fTTI i, the 
fermionic system operator d(t) e \a(t), a^(t)} at time t only 
depends on the input noise operators at time t' < t. As a re- 
sult, in the Markovian limit, the causality relation read s 30 ' 32 

[Xin if') ,d(t)] + = 0, for t' > t. (24) 

For the similar reason, the system operator at t is independent 
of the output noise operators at time f < t 

[£ou, (f'M(0] + = 0, forr'<r. (25) 

According to Eqs. (l22H25b . the anti-commutators between 
noise and system operator are converted to those between sys- 
tem operator a 30 ' 32 

[£ in (t'),d(t)l = y L 6(t-t')[a(t'),d(t)] + , (26a) 
[ft in (f'M(0] + = y R 0(t-t')[a(t'),d(t)] + , (26b) 

where step function 6{f) is defined as 

(1, t >0 
6(t) = \ i, t=0 . (27) 

(6, t < o 



For the multi-level case, this causality relation Eq. d26l i can be 
generalized to the system projection operators cry, i.e. 

[Xin(0>^7«] ± = 7L0(t-t')[a(t'),cnj(t)] ± , (28a) 
[ft in (f'),^y( f )] ± = 7R0(t-t')[a(t'),o- u (t)] ± . (28b) 

The choice of the commutative and anti-commutative relation 
in Eq. ( l28l is determined by the parity of the electron number 
change, see Eq. ( TT3T >. 

In the following, to simplify the notation, we will omit the 
subscript "in" of the input noise operators, unless stated oth- 
erwise. 

D. Current and current noise 

For the quantum transport problem, we are interested in the 
average current and the current noise spectra. In this subsec- 
tion, we will give the expressions of such quantities in terms 
of the noise operators. 

We consider the current through the right lead for example. 
For simplicity, let us first study the single level case. The 
formula for the multi-level case with Coulomb blockade will 
be discussed later. The current operator can be defined as the 
changing rate of the electron number on the right lead, i.e. 

Ir = jNR = Yu£j c ) a + h - c - (29a) 

j 

= y R a* (t) a (t) - ft 1 (t) a (t) - a 1 " (f) K (t) . (29b) 

The second line is obtained by noticing the relations in Eq. © 
and Eq. ((9). We point out that the current operator can be di- 
vided into two parts: (i) the damping part jro^ (t) a (f), which 
is proportional to the level occupation and the escaping rate 
jr, and (ii) the fluctuation part (the last two terms), which is 
induced by the noise operators H (t) and fC (t). 

For the average current, we take the average of the current 
operator Ir over the thermal states of the leads 

(t R ) = y R {a? (t) a (f)> - <ft f (0 a (f)> - (a 1 " (f) % (f)>. (30) 

And for the current noise, we first calculate the current-current 
correlation function 

g {1) (r) = lim Re \{I R (t) t R (t + r))l - (I R f (31) 

At steady state, its Fourier transformation gives the current 
noise spectrum^ 

Jr*co 
g (2) ( T ) cos (cot) dr. (32) 
o 

To calculate the correlation (I R (f) I R (f + t)) in g (2 ' (t), by the 
definition of I R in Eq. d29bl ). one need to calculate the two 
time correlation such as 

(a + (f) a (t) a 1 " (t + t) a (t + t)>, (33a) 
(a* (t) <R it) a 1 " (t + t) a (t + t)>, (33b) 
(a 1 " (f) n {t) 1^ (t + T)a(t + t)>. (33c) 
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In Sect. [En] we will show that the fluctuation part in the cur- 
rent operator does not contribute to the average current, so the 
average current (/«) = y«(a' (f) a (f)) is held. But the fluc- 
tuation terms will contribute to the current noise through the 
correlations in Eq. d33b . 



III. APPLICATION I: SINGLE LEVEL TRANSPORT 

In this section, the general quantum Langevin formula is 
applied to the resonant transport through a quantum dot. As 
the first example, we consider a model in which only one sin- 
gle energy level in the quantum dot is relevant. The system 
Hamiltonian reads 



tiojaa' a. 



(34) 



We consider the large bias condition and assume that the 
single-particle energy level with energy ha>o is well within the 
bias window, i.e. fii - cl>o, wo - fiR » yi, jr, for be- 
ing the chemical potentials of the left/right leads. According 
to the discussion in Sec. Ill Al the quantum Langevin equation 
reads 



5(f) = - jL + yR a(t) + Z(t)+ft{i), 



where a it ) = a (t) e ilJ <", £ (t) = e ioJo '£ (f), and "k (t) = e icJo "R (t) 
are defined in the rotating reference frame to single out the 
slow-varying dynamics. In the white-noise limit, the correla- 
tion between the noise operators can be written as (see Ap- 
pendixlAli 

UHO 1(f)) = 7 L Sit-lf), (36a) 
<*(/)#(0> = y R 6{t-f), (36b) 
(£{t)l'{t')) = <ft f (*)&(*')> = 0. (36c) 

Using these relations, we calculate average current and current 
noise. 



A. Average current 

From Eq. ( 1351 . the system operator a (f) in terms of the 
noise operators is 

a(t) = e- r -'a(0)+ f e-^'-'^KOdt' 
Jo 

+ [' e-^'- t ' )i R{t')dt', (37) 
Jo 

where T - Jl + Jr- Multiplying the noise operator £? (t) on 
both sides of Eq. (fJTJ, we have 

{ZHt)a(t)) = e- T -'UUt)a(0)) 



+ [' e-^-''\lHt) l{t'))dt' 
Jo 



n 

2 ' 



(38) 



Here we have assumed that at initial time t — 0, the sys- 
tem and the reservoir are independent, i.e., {JJ (t) a (0)) = 
(j? (t))(a (0)) = 0. Similarly, we obtain 



(a t a)X«> = y, 



and 



<5 f (f)9?(0> = <ft T (f)a(r)> = 0. 



(39) 



(40) 



Thus, according to Eq. ( |29bb , the fluctuation part of the cur- 
rent operator does not contribute to the average current, and 
the average current becomes 



(Ir) = y R (a a). 



(41) 



In order to determine the mean occupation number {a 1 a), we 
use the equation of motion 



dt 



-a'd — d^a + a^a 



= -ra ' a + a T (t) z (*) + £ (?) a (t) 

+aHt)K{t)+ft' {t) a {(). (42) 



(35) 

The ensemble average leads to 



d , 

—(a' a) = -T{a'a) + y L . 
dt 

Thus, the averaged population in the quantum dot is 



(43) 



<atg> = — y ± Z^_ e -(ri+r«)'. (44) 

yi + r« yl + yr 

As a result, the average current at steady state for t — > +oo is 

YlYr 



(h)s 



Yl + Yr 



(45) 



which is the well-known result for the resonant tunneling 
transport Aii 



B. Current noise 

To investigate the current noise, we calculate the current- 
current correlation (Ir (t) Ir (t + r)). With the definition of the 
current operator in Eq. d29bt , the noise contains typically two- 
time correlations like 

(a 1 " (f) a (t) a 1 " (/ + t) a (t + r)) = (ft (t) h(t + r)>, (46) 

and 



<a T (f)R(f)R T (f + T)a(f + T)>. 



(47) 



We will discuss such correlations one by one. 

Noticing that the electron number correlation function 
(h (t) h{t + t)) contains only the system operators, we use the 
quantum regression theorems! and Eq. d43l and obtain 



— (h (t) h(t + r)> = -T(h (t) h(t + T )> + y L (h (0>. (48) 
dr 
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This equation, together with the initial condition with respect 
to r, i.e. for r = 0, (ft (t) hit + t)> = (ft (t) h (/)) = (ft (£)), 
determines the occupation number fluctuation in the quantum 
dot. The steady state correlation is 



lim (h (t)h (t + t)> 



JLJR 



-(7l+7r)t 



(7l + Yr) 2 ijL + 7r) 2 



(49) 



The other terms contain the correlations between the system 
and noise operators. Taking (5 1 (t) % (t) W (t + r) a (t + r)) for 
example, according to Eq. Q71 l. we have 

(a' (t)R(f) $ (f + t) a (t + T)) 

r-t+T 

dt x dt2e-^'-' ,) -^' +T -' 2) G(t l ,t,t + T,t 2 )A50) 



Jo Jo 

where the four-time noise correlation is defined as 

G(/i , h, f 3> f 4 ) = <X f (h)K (t 2 ) ® (t 3 ) I (t A ) >. (51) 

According to the independent noise assumption and the white- 
noise approximation, 

G(t u t 2 , h, f 4 ) = y L y R 6 (h - f 4 ) 6 (t 2 - h) . (52) 

Thus, we have 

{aHt)ft(t)ftHt + r)a(t + t)> 

Xt r-t+T 
dh \ dt 2 [e-^ 2,+T - h - h) 6(h -t 2 )5(T)\ 

= e 2 o(t), for/— > +oo. 

JL+7R 



(53) 



Similarly, 



(aJ (/)??(/) a f (/ + T) a (f + t)> 



< 



-(n+y s )T for f +c 



(54) 



It can be checked that all the other terms in the current-current 
correlation function vanish. Consequently, the current-current 
correlation function is 



(r) = - 



2 2 

r L r R 



_ e -rr + _Xj2*_ g -ST tf(T)> (55) 

(jL + JrY JL + JR 



and its Fourier transformation gives the current noise spectra 

rl + y 2 R + " 2 



S (to) = 2e(I R ) s 



(jL + JR) 2 + 



(56) 



This result accords with the ones derived from other 
approaches, 17 and shows that the presence of the single level 
quantum dot suppresses the zero-frequency current noise to 
half of the Poisson value S p - 2e(/«),j in the case yi = jr. 

It is worth to emphasize that, clearly shown in our quantum 
noise approach, although the fluctuation part [see Eq. d29bH 



of the current operator does not contribute to the average cur- 
rent, it does the current noise. According to our approach, 
the current-current correlation originates from three different 
kinds of sources: (i) the on-site number-number correlation 
[Eq. j49ll, which always contributes a positive correlation, (ii) 
the correlation between the fluctuation terms Eq. d53"l l, which 
induces a white-noise correlation, and (iii) and the correlation 
between the on-site number and the fluctuation term Eq. d54b . 
which always provides a negative correlation. This classifica- 
tion of current-current correlation is also valid in the interact- 
ing case, as will be discussed below. 



IV. APPLICATION II: COULOMB BLOCKADE 

A. Average current 

Now we apply the general theory to the Coulomb blockade 
case. For simplicity, we assume that only one single orbital 
level in the quantum dot is relevant (i.e., within the energy 
range of interest). The system Hamiltonian reads 

#sys {o-i, — ftu>ia^ai + Tzw^ajai + t/a|a-|-a|fl|, (57) 

where Tilj-\^ are the single electron energy for spin-up and 
spin-down electrons in the quantum dot, and U is the Coulomb 
interaction strength between two electrons. In this paper, we 
consider the large U limit, i.e. tiio-^ + U, Ha>i + {/»//£» 

As has been discussed in Sec III Bl though there is interac- 
tion between the electrons in the system Hamiltonian, H sys is 
diagonalized as 

H sys = ftco^cr f-f + fi«Lcr u + (ftw T + hu l + U) cr dd , (58) 

and the projection operators are related to the Fermion opera- 
tors by 



a[ = O-yJ, + a- rd , 



(59a) 
(59b) 



where cr l; - = for i, j = v, 1\ I and d. The subscripts 
v, t , I and d represent the vacuum state, spin-up, spin-down, 
and doubly occupied state, respectively (Fig. [2j». Annihilating 
an electron with definite spin (say spin-up) from the quantum 
dot consists of two different projection processes depending 
on whether the spin-down level is occupied or not. 

Here, we assume that the quantum dot is coupled to fer- 
romagnetic leads. Thus, the electron with different spin can 
tunnel on and off the quantum dot with different rates. The 
quantum Langevin equations of the projection operators o",y 
in this Coulomb blockade case follow the general formula in 
Sec. HI Bl The resultant equations for the diagonal elements 
are 
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and those for the off-diagonal elements are 



dvr = 


-j'A vT cr vT 


0"J.<i = 


-iA id cr id 


&vl = 


-zA v jo- v j 


o~1d = 




& vd = 


-iA vd cr vd 


& n = 


-iAij.tr |j 





+ Ti 




1 

2 


r t 


+ T| 




2 


r T 


+ Tj 




1 

2 


r T 


+ Tj 




2 


r T 


+ Tj 




2 


r T 


+ Tj 


2 



— tr vT - [T T (trjj + tr^) - Tjtr iT - T^o-^/] , 
-°"vj. - yO"Trf + [^i (°"vv + o"u) + ^T°"n - yffvd] » 



T I Tt r +"i 

— ~ cr ]d + —tr.j + [:Fj (tr TT + cr dd ) - ^ T tr n + tr vd J 



(60a) 
(60b) 
(60c) 
(60d) 

(61a) 
(61b) 
(61c) 
(61d) 
(61e) 
(61f) 



r 



where Ay = ojj - w,-, f s = T s (t) = £ s (t) + %,(t) for s e {|, 1), 
and the spin dependent noise operators -C s {t) and 7?j(f) are 
defined as 

£ s (t) = -J^-^'MO), (62a) 

ft 

«,(/) = -^c,t' /,,V '<>(()). (62b) 

The damping rate r. s = ji s + y« s , with the spin dependent 
tunneling rates 



y Ls = 2nf s (a> s )Df ) (aj s ), 
y Rs = 2ng(u> s )Df ) (a> s ), 



(63a) 
(63b) 



where f s and tf s are the coupling amplitudes of the quantum 
dot to the left and right leads, and D^\a>) and (cS) are the 
spin-resolved density of states of left and right leads, respec- 
tively. 

These Langevin equations of the system variables are anal- 
ogous to the ones used to describe the quantum theory of 
Lase r 23 ' 24 . In the quantum theory of Laser, the atoms are sub- 
ject to bosonic reservoirs, while in our quantum transport case, 
the quantum dot is "pumped" by a fermionic reservoir (the 
left lead), and output to another fermionic reservoir (the right 
lead). 

In contrast to the non-interacting case [see Eq. d35l)l, the 
noise operators couple to the system projection operators in 
Eqs. ( l60l l and ( f6Tb . The correlations between noise operators 
and projection operators, such as (-Ci(t)<Tij(f)}, are calculated 



according to the generalized causality relation Eq. d28| i. Tak- 
ing (X|(f)o"vT( f )) for example, 



<4«<M(0> = (X|(f)tr, T (0) 

-y L ([a\(t), ovtCOU) = \-y L <tr n .(f) + tr TT (f)> , (64) 



where tr v ,j(f) = tr,,j(f)e" 1,T ' is the slow-varying amplitude of 
projection operator, and Xj(f) = X|(f) e '" T '- The correlation 
(tr v -|-(f)-£|(f)) m me second line of Eq. d64l i vanishes, when the 
noise operator Xj(f) acts on the full-filled fermi sea of the left 
lead, see Appendix lAl 

With these correlations, ensemble average of Eqs. d60l l and 
(l6"TT l gives the "rate equations" for the diagonal elements: 



<o"w> = - (tlt + 7u) <o"w) + TrtW) + r«i<o"li)(65a) 

(d- TT ) = -y RT <tr TT > + jL](o- V y) + (y L { + Jr{) (a- dd ),{65h) 

= -7/u<0"u> + ru<°"iv) + (tlt + Trt) <o"^>,(65c) 

(d-dd) = - (yz.T + Jn + 7«T + YRl) (°~dd), (65d) 
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and for off-diagonal elements: 



(0"vt) = 


- 0'A„t 


+ ri)<cr l , T ) - 


7U + 2jRi , . 

— ~ — <0V>, ( 66a ) 

2 


/ * \ 




+ 1 2) {a- id) ~ 


— <cr„ T ) , (66b) 


(irvO = 


-(iKi 


+ r 3 ) (o-vi) + 


2 (°W> ( 66c ) 


(°~fd) = 


- (i^d 


+ r 4 ) (o- v ) + 


yW>. (66d) 


(cr v d) = 


- (iAvd 


+ 2y LT + 2y u 


+ 7rt + 7*1 ) (°"vd> .(66e) 


<o"Tl> = 


- ('A T j, 


+ JR] + Jn) < 


o- n ), (66f) 



where Y x = (y LT + y u + y RT )/2, T 2 = {y L ] + y L \, + y«T + 
2y Rl )/2, T 3 = (y LT + y u + y R |)/2 and T 4 = (y LT + y u + 2y R] + 
yui)/2. Eq. d66l l shows that the coherence between energy 
levels vanish after a long time, i.e. 



((Tij(t)) = 0, for i + j, and t — > +oo. 



(67) 



This indicates the two different spin channels are incoherent, 
which physically arises from the fact that noise operators with 
different spins are uncorrelated, i.e. (Xl(0-Cj.(f')) = 0. 

The rate equations d65l l describe the population transfer be- 
tween each energy levels, and the steady state populations are 



(cr n ) 
<0"U> 

(o-dd) = 
The average current is 



7r]7ri 


jL{jR-\ 


+ 7l-\7ri 


+ 


7RTTRI 




Jl\7ri 








+ 7L17RI 


+ 


7R17RI 




7li7ri 






7L17R1 


+ 7LT7RI 


+ 


7RT7K1 



(68a) 

(68b) 

(68c) 
(68d) 



Or) = </« T > + ihi) 
y ST (cr TT ) + y«x<cr u ) 

(yiA+7Li)7RWRl 
yuy«T + TiTTRI + 7R17RI 



(69) 



The current vanishes if y«| = or y^j, = 0. This is because 
turning off a certain spin channel, say the spin-up channel, i.e. 
ysj = 0, will induce the accumulation of the spin-up elec- 
tron on the quantum dot. Then the electron tunneling of both 
spin channels is blockade due to the strong Coulomb inter- 
action. When the tunneling rates are spin-independent, i.e. 
7li = 7U = 7l and^y ST = y Rl = y R , the average current in 
Eq. d69l ) becomes (/«) = 2yiy R l{2yi + y«), which accords 
with the results obtained by other methods i 11 - 16 ' 17 



B. Current noise 

Now we turn to the current noise. Similar to the non- 
interacting case, the two-time correlations of the following 




FIG. 2: (color online) Schematic illustration of Coulomb blockade 
transport. The quantum dot can be regarded as a four-level artificial 
atom with |v), | t)> I -I) an d \d) representing the vacuum, spin-up, 
spin-down, and the doubly occupied states. The lower panel shows 
the two element processes of a spin-up electron tunneling out of the 
quantum dot depending on whether the spin-down state is occupied. 



form should be calculated 



(n s (f) n s , (t + t)>, 

<aj(0 + t) a, (f + r)), 

{a\ (r) n s (t) a\, (t + r) a s , (t + r)>. 



(70a) 
(70b) 

(70c) 



In the interacting case, the system projection operators can- 
not be expressed in terms of the simple integration of the 
noise operators as in the non-interacting case, see Eq. d37l ). 
The causality relations introduced in Sec III CI provide us a 
convenient way to convert the noise-system correlation to the 
system-system correlation. Thus, in the white-noise limit, as 
a powerful tool, the quantum regression theorem is applied to 
calculate the two-time system correlations. 

Noticing that the noise operator 'Rs(t) plays the role 
of "annihilation operator", the correlations between noise 
and system operators can be calculated following the 
spirit of the Wick's theorem(see Appendix [A}. Tak- 
ing the spin-up component for example, the correlation 
(a\ {t) 7? T (t) n\ (t + t) a T (t + t)> is 

(a] (t)K T (r)*J (f + r)a T (f + t)> 
= {a\(t)VR^t),n\{t + r)] + a^t + T)) 

-{a\ {t)<R\ (t + T)R t (t)a t (f + r)) 
= y R (a\(t)a r (t + T))6(T). (71) 

The second line of Eq. < T7TT > is simplified by noticing the fact 
[7?t (t) , ?£| (f')]+ = 7rSU - 1'), and the third line vanishes since 

[a| (f) , ??| (f + t)]+ = 0. This white-noise correlation provides 
a constant current noise background. Due to the 6(t) function, 
only equal-time correlation (r = 0) is relevant. By noticing 
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Eq. ( 1691 , this correlation is written as 

(a\ (f)K T (t)R] (t + T)a T (t + r)> = </ rt )(5(t). 



(72) 



For the correlation (aj (f) ^ T (f) aj (f + r) a T (r + r)> = 

<a| (O^t ( f ) n T ( f + T ))' ^ can be translated into the correla- 
tions between the system operators using the causality rela- 
tions 

{a\ (0« T W»T ( f + ^)> 
= {a\(t)VR^t)^(t + T)]) 
= y R {a\(t)[a^{t),n^{t + T)}) 

= y«<n T (f) « T (r + r)> - y R {a\ (f) « T (r + r) a T (f)>. (73) 

The first term cancels out the contribution of Eq. ( |70at . As a 
result, the spin-up current-current correlation is 

*® (t) = Hm </ ST (0 /jet ( f + T )> 
= ; lim(/ ST )5(T)+r2 T < fl i( f ) nT ( f + r)a T (0). (74) 

The current-current correlations of different spin components 
are calculated similarly, and in general, they can expressed in 
terms of the correlations of system operators as 

gfl{T)= )im (l Rs (t)l R At + t)) 

r-»+oo 

= lim (I Rs )6(t)8 ss > + y Rs y Rs ' (a\ (t) n s > (t + r) a s (f)>(75) 



<D 



=3 

o 
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Here we have shown an analogous form of the current-current 
correlation to the second order optical coherence function, 2 - ; 
The last term of Eq. d75l l can be calculated from the quantum 
regression theorem. By this theorem, the current-current cor- 
relation function is determined by the rate equations (l65l ). and 
in the Coulomb blockade case, it does not show the effect of 
the quantum coherence terms in Eq. d66i l. 
The total current correlation function is 

8 {2) (t)= lim(i R (t)I R (t + T))-(I R ) 2 

t— >+oo 

- ~ (2) (t) + gfl (t) + gf] (t) + gfl (r) - (I R f (76) 



'T1 



'IT 



'11 



Its Fourier transformation gives the current noise spectrum 
S(a>). In the spin independent tunneling rate case, i.e. y^ - 
la = 1l and y ST = y Rl = y R , the noise spectrum is 



S(d) = 2e(I s ) 



4yj + 3y L y R + y\ 
(2y L + y R ) 2 + or 



(77) 



This result deviates from the single-level case [see Eq. (|56l)l, 
due to the presence of the Coulomb interaction. 

Typical current noise spectra for the spin dependent tunnel- 
ing rate case are shown in Fig. ([3^). The Fano factor is 

F = S(u> - 0)/2e(I R ) 

j _ ysty«i(ygy«T + t^trQ - yma(y^ - ^i) 2 ^ 
2(yLiy«t + ymn + yfity«i) 2 



FIG. 3: (color online) (a) Current noise spectra S (a>) (normalized by 
the Poisson value S P = 2e{I R )) for y Rl = 0.1,0.3, • • ■ ,0.9. Other 
parameters are chosen as y^ = y L ± = y^ = 1. (b) Fano factor as 
a function of the imbalance between spin-resolved tunneling rates 
P L and P R , which are defined as P L = yi^KyiA + 7u) an d Pr - 

y«T/(y«T + y«i)' for g iventotaltunnelm g rates 7'iT+ru =y«T+rn = 

1. The white thick lines are the boundary between sub-Poisson and 
super-Poisson regime, i.e. F = 1. 



It is found that super-Poissonian noise arises when the tun- 
neling is spin dependent, which can be realized, e.g., by us- 
ing magnetized barriers between the leads and the quantum 
dot. Super-Poissonian noise appears when the numerator of 
the second term becomes negative. The Fano factor as a func- 
tion of the tunneling rate imbalance is shown in Fig. ^p). 

Physically, the super-Poissonian noise is the consequence 
of the dynamical channel blockade effecti^^ The tunnel- 
ing rate imbalance induces different average currents for the 
two spin channels. Thus, in additional to the noises of each 
channels themselves, the shot noise between the two channels 
gives rise to the low frequency noise enhancement. Such kind 
of shot noise is absent when Pi — 1 - P R since the two spin 
channels have the same current. 
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V. RELATION TO OTHER THEORIES 
A. Relation to Landauer-Biittiker formula 

Here, we show that the Landauer-Biittiker formula can be 
reproduced by the quantum Langevin approach. For simplic- 
ity, let us consider the single energy level transport example. 

According to Eq. (|29bt and the boundary relations Eq. d23l ), 
the current operators can be expressed solely by the input and 
output noise operators. For example, 

Ir = — (Kt W ^ut (0 - $1 it) Kn (/)) . (79) 

Thus, it is clear that the average current is divided into the in- 
put current proportional to 0tL (t) Ki n (t)) and the output cur- 
rent proportional to 0V oyA (/) ^ ollt (t)). 

Furthermore, defining the scattering matrix S, the Fourier 
transformation of output noise operators is expressed in terms 
of the input operators as 



Xout (to) 
Kut (to) 



= S 



■Cm (to) 
#in (W) 



with 



-£in/out (to) 

^in/out (to) 



e iu 'l in/out (t)dt, 

CO 



CX 



(80) 

(81a) 
(81b) 



and 



S(a>) = 



T r ^ l (lo) R R ^ R (to) 

JL-JR 



Jl+Jr- 2i(0 



JR 



7l 

Jl-Jr 
2 



+ ILO 



,(82) 



where the functions T^j(to) and R,^ ; (<y) can be regarded as 
the energy dependent transmission and reflection coefficients 
from lead / to lead ;. The Fourier transformation of the aver- 
age current is 

(t R (to)) = J" {h(t))e ia 'dt 

i r°° 

= — [<^ ut K)Kut(w' + w)) 

d( / 

- (Kl(a>')K m (a>' + (o))]— . (83) 

Noticing the relation Eq. ( f80b and the correlations between the 
noise operators 

Ul(to')£ m (to' +co)) = 2ny L 6(co), (84a) 
(Kl(co')ft m (co' + co)) = 0, (84b) 

we obtain the Landauer-Biittiker-like formula of the average 
current 



(I 



r) = T{co )— 

J-oo 2n 



IlJr 

JL+JR 



(85) 




20 30 
Frequency m 

FIG. 4: (color online) Local density of states (LDOS) for the 
Coulomb blockade case obtained by the quantum Langevin ap- 
proach. Parameters are a>f = 10, V = 30 and y = 1. LDOS shows a 
double-peak structure for U » y. 



with the transmission spectrum 
T(co) 



JhjR 



(ZilZ«)2 + w 2- 



(86) 



B. Relation to non-equilibrium Green's function theory 

Here we discuss the relation between the quantum 
Langevin approach and the non-equilibrium Green's func- 
tion theory for the quantum transport problems. The retarded 
Green's function is defined as2i 



G s (T) = -iG(T){{a s (t + T),al(t)}), 



(87) 



for s =| or X, from which the local density of states (LDOS) 
D s (a>) is given by 



D s (co) = ~~hn[G s (co)l 
n 



(88) 



where G s (a>) is the Fourier transformation of G s (t). The 
LDOS contains the essential information about the system rel- 
evant in quantum transport. In the following, we take the 
Coulomb blockade example, and give the retarded Green's 
function and the LDOS using the quantum noise approach. 

Noticing that the definition of the retarded Green's function 
Eq. d87l i only involves the system operators a s (t) and a\(0), 
we apply the quantum regression theorem to calculate their 
correlations. The retarded Green's function can be expressed 
in terms of the two-time correlations between the projection 
operators. Consider the spin-up component for example, 

-i9(f) <{cr„ T (f + t) - cr ld (t + r), 
o>(0 - o- dl (t)}) . 



G T (T) 



(89) 



The equations of motion for these projection operators are 
given in Eq. d66l i. By the quantum regression theorem, the 
two-time correlations are determined by 



d_ I (o- Tv (i)o- vt (t + t)> 
dT \ (o-t v (t)cr ld (t + t)> 



M 



<cr Tv (f)o- v , T (t + t)> 
(cr rv (t)o- ld (t + r)> 



(90) 
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with the initial condition for r = 



(o>(*)o- vt (t)> \ _ / <o- TT (r)> \ 
(^ v (t)cr ld (t))J \ j' 



where the coefficient matrix M is defined as 



M = 



-3y/2 - ;<j>t 



-3y/2 



-y/2 -57/2 - z(to T + 17) 



(91) 



(92) 



Here, y^ = y^ = 7^ = 7^ = 7 is assumed for simplicity. 
The other correlations involved in Eq. ( 1891 can be similarly 
calculated. 

Thus the retarded Green's function is 



(93) 



G t (t) = -iB{T)e- lyT (w + e-^ r + W-e- ia, - T ) , 
with the renormalized frequencies 

u ± = w T + ~ ± i ^U 2 -2iyU-4y 2 , (94) 

and the weight factors 
1 



W ± = - 



C//6 - iy 



2 2 Vt/ 2 - 2/7!/ - 4y 2 



(95) 



The Fourier transformation of the Green's function gives the 
LDOS (see Fig. [4). It is obvious that, for the large U case 
considered in this paper, the LDOS consists of two Lorentz 
shape peaks, centered around aif and a>f + U. The two peaks 
separate from each other by U, which is a signature of the 
Coulomb blockaded 



VI. CONCLUSIONS AND OUTLOOKS 

In this paper, we have developed a quantum noise approach 
to treat the quantum transport through a nanostructure such as 
a quantum dot. We formulate the average current and the cur- 
rent noise in terms of the correlations between the noise opera- 
tors. The quantum noise approach is applied to a paradigmatic 
example, namely, transport through a single quantum dot un- 
der large biases and both the non-interacting and Coulomb 
blockade cases are investigated. With the Markovian approx- 
imation for the tunneling processes, the electron-electron in- 
teraction in the quantum dot can be exactly treated. 

The quantum noise approach provides a bridge between 
quantum optics and quantum transport. Thus notions and 
methods in the quantum optics could be adopted to study 
quantum transport through nanostructures. Although we show 
the application of the quantum noise approach by a single 
quantum dot example, the theory is not limited to this simple 
case. On one hand, the system could be generalized to more 
complicated ones, such as coupled quantum dots, multi-end 
nano-circuits, or systems with spin interaction. On the other 
hand, the reservoirs of other kinds, such as phonon baths or 
spin baths, could be included to explore how such reservoirs 
would affect the current and current noise, providing a method 
of studying the bath dynamics via current noises. The Marko- 
vian approximation may also be released with colored noise 
correlation functions of the reservoir used in lieu of the white- 
noise model adopted in this paper. 



APPENDIX A: PROPERTIES OF THE NOISE OPERATORS 

In this appendix, we give the correlations between noise op- 
erators. We consider the single-level case here. The physical 
quantities of interest are determined by the noise correlations 
such as (£"'" (f) -C (?'))• According to the definition of the noise 
operators, 

<r(ox(o> 



= 2 te^-^'-^-^ {b\b k .) 

k,k' 



f(co k )D(co k )nf (fi> k )e«*- aM ~ f >d<o k , (Al) 



where D(a>k) is the density of states in the leads, and 

1 



l _|_ e (tiot-fi L )/k B T ' 



(A2) 



is the thermal occupation number of the lead in quasi- 
equilibrium. The Markovian approximation requires two as- 
sumptions. First assumed is the "flat band" condition that the 
relative change of the effective density of states around the 
resonant a>o over a range of the characteristic damping rate jl 
is much less than unity, i.e., 



d In D(cj k ) Y 
du) k J 



» JL, 



(A3) 



where D{<jj k ) = g 2 (a) k )D(a) k ), Under this condition, D{u)k) can 
be replaced by its value at ojq, and the correlation becomes 



f 

Jo 



Wt-GJoX'-O 



da>k- 



(A4) 



Here, the zero temperature case has been considered for sim- 
plicity. Second, under the large bias condition, the resonant 
level a>o is far away from the fermi energy and the conduction 
band bottom (chosen as the energy origin), i.e. 



li L - co , oJo » 7l- 



(A5) 



In this case, the integration over a> k is extended to ±00, and 
finally results in the white-noise correlation 



UHt)£(t'))=y L 6(t-t'), 



where y L = 2n^ 2 (a> ) D (cj ). 
Similarly, for the right lead, 



<£(t)#(0> = yji«y(f-0. 



(A6) 



(A7) 



Here, we have use the fact that the thermal occupation number 

n ttP {^j) ~ f° r riS^t l ea d ^ound the resonant level ojq. 
In the same way, one can show that other noise correlations 
vanish, i.e. 



<£(*)£ + (O> = <£ + «£(O> = 0. 



(A8) 
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Note that Eq. ( IA8b implies that the noise operators X) (t) 
and % (f) play the role of "annihilation operators", since they 
always give zero correlations when they stand on the right- 
most position. With this observation, the normal-ordered 
product of noise operators can be defined by placing £} (f) 
and ft. (f) to the rightmost position, and the expectation value 
of the normal-ordered product vanishes identically. Thus, 
the Wick's theorem is generalized to the noise operators and 
the current and current noise can be exactly calculated in the 
white-noise limit. 
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